Disorder-induced superconductivity in ropes of carbon nanotubes 
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We study the interplay between disorder and superconductivity in a rope of metallic carbon 
nanotubes. Based on the time dependent Ginzburg Landau theory, we derive the superconducting 
transition temperature T c taking into account the critical superconducting fluctuations which are 
expected to be substantially strong in such low dimensional systems. Our results indicate that, 
contrary to what is expected, T c increases by increasing the amount of disorder. We argue that 
this behavior is due to the dynamics of the tubes which reduces the drastic effect of the local 
disorder on superconductivity by enhancing the intertube Josephson tunneling. We also found that 
T c is enhanced as the effective dimensionality of the rope increases by increasing the number N of 
the tubes forming the rope. However, T c tends to saturate for large values of N, expressing the 
establishment of a bulk three dimensional (3D) superconducting order. 



INTRODUCTION 

Since their discovery in 1991 [1], carbon nanotubes 
(CNT) have been studied under close scrutiny due to 
their eye-catching properties which are of a great inter- 
est not only for nanotcchnology but also for fundamen- 
tal physics. A carbon nanotube, which can be regarded 
as a tiny cylinder rolled up from a graphene sheet, is 
a good candidate to study electronic properties in one- 
dimensional (ID) systems where electron-electron inter- 
actions are substantially important. 
CNT can be synthesized as a single walled tube (SWNT) 
or multiwallcd tubes (MWNT) consisting of two or more 
concentric shells. SWNT can also be assembled to form 
ropes of ordered parallel tubes arranged in a triangular 
lattice [2-4]. The nanotubes, which are nearly of the 
same diameter, can have different kind of hclicitics, but 
in general | of them are metallic [5, 6]. 
The transport properties of the rope is found to be 
strongly dependent on the amount of disorder within the 
tubes [7, 8] . It has been reported that the intertube elec- 
tronic transfer is enhanced in the presence of disorder, 
leading to a charge carrier derealization [9] . This feature 
raises the question whether such disorder-induced inter- 
tube coupling can be observed for the superconducting 
order in ropes of CNT? 

The first superconducting signature was observed, in 
1998, as a proximity effect in isolated metallic bun- 
dled SWNT connected to superconducting leads [10, 11]. 
Later on, intrinsic superconductivity has been reported 
in ropes of CNT with a transition temperature T c = 0.55 
K [6, 12, 13]. Ferrier et al. [6] studied the dependence of 
the superconducting transition temperature on the num- 
ber of the metallic tubes included in the rope and on the 
amount of disorder. They found that superconductivity 
arises only in ropes with more than 100 tubes. Bowever, 
the most striking result of Rcf.[6] is that disorder, con- 
trary to what is expected, may induce superconductivity: 
the larger the amount of disorder, the stronger the su- 



perconducting correlations. Nevertheless, at a very large 
disorder amplitude, the superconducting order collapses 
as in other superconducting materials. 

Superconductivity at T c = 15 K has been also reported 
in zeolite-inserted SWNT of small diameter (0.4 nm) [14]. 
Takisue et al. [15] found a superconducting transition at 
T c ~ 12 K in MWNT encapsulated in zeolites. These 
relatively high critical temperatures put the question on 
the origin of superconductivity in SWNT. Bow can a su- 
perconducting order develop in such low dimensional sys- 
tems where thermal fluctuations are expected to destroy 
any long range ordered state? The surprising observation 
of superconductivity in CNT has stimulated many theo- 
retical studies to found out the underlying mechanism. 

The realization of a superconducting order in ropes of 
CNT has been ascribed by Gonzalez [16] to the presence 
of strong attractive electron-electron interactions medi- 
ated by phonon exchange. The latter prevails over repul- 
sive Coulomb interaction in ropes with hundred or more 
of metallic nanotubes. 

Other models based on phonon mediated attractive 
mechanisms have been also proposed [5, 17, 18]. In par- 
ticular, the dependence of the superconducting transition 
temperature on the number of tubes was quite under- 
stood in the framework of the model elaborated by Eg- 
ger and De Martino [5, 18] who introduced the Josephson 
couplings between the tubes and the phase fluctuations 
of the superconducting order parameter. Bowever, a pro- 
nounced discrepancy with the experimental data emerges 
with decreasing the number of the tubes embedded in the 
rope [19, 20]. 

To explain the relatively high superconducting criti- 
cal temperature reported in SWNT, Sasaki et al. [21] 
have proposed a new mechanism where superconductiv- 
ity originates from the edge states specific to graphene. 
The authors argued that superconductivity is due to a 
superconductor/normal/superconductor junction where 
the superconducting phase is realized at the ends of the 
SWNT while the bulk part of the tube remains metallic. 
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An other scenario has been proposed by Zhang et al. 
[22] to account for the occurrence of superconductivity 
in SWNT connected to superconducting or normal elec- 
trodes. The authors argued that the SWNT becomes 
superconducting in the range of 11-30K due to the pres- 
ence of van Hove singularities in the electron density of 
states of the nanotube. 

Karnaukhov and Diks [23] ruled out the electron- 
phonon interaction mechanism to explain the formation 
of the superconducting state in SWNT due to the rela- 
tively large value of the critical temperature. The authors 
suggested an alternative attractive electron-electron in- 
teraction originating from strong hybridized interaction 
induced by the two-band electron structure of SWNT. 

Recently, Belluci et al. [24] have theoretically ar- 
gued that superconductivity can arise by a purely elec- 
tronic mechanism in ultrasmall diameter SWNT and end- 
bonded multiwalled ones due to the screening of the for- 
ward scattering processes. 

More recently, Le Hur et al. [25] have derived a theo- 
retical model to study the possibility of a superconduct- 
ing proximity effect in metallic SWNT in the presence of 
superconducting substrate. The authors showed that the 
latter induce an unconventional double superconducting 
gap in the tube. 

The outcome of the above-mentioned studies is that 
the origin of superconductivity in CNT based systems is 
still under debate and many relative issues are not yet 
totally unveiled. In particular, the role of disorder on 
the stability of the superconducting phase has not been 
addressed in previous theoretical studies [20]. This is a 
key point which may shed light on the formation of the 
superconducting phase in low dimensional systems. 

In this paper, we theoretically investigate the effect of 
disorder on the superconducting state in a rope of CNT. 
The model is based on the time dependent Ginzburg- 
Landau theory taking into account the superconducting 
fluctuations which are substantially important in CNT 
regarding their low dimensionality Ferrier et al. [19, 26] 
have actually observed, in ropes of CNT, a large do- 
main of superconducting fluctuations which extends to 
IK, namely twice the transition temperature (T c =0.5 K). 
In the following we present our model and discuss the ob- 
tained results in section III. Section IV is devoted to the 
concluding remarks. 



simplicity we assume that all the tubes are metallic while 
experimentally |, on average, are semiconductors. This 
assumption does not affect the outcomes of the present 
model which depends basically on the amount of dis- 
order in the rope and on the intertube Josephson cou- 
plings. From the numerical point of view, one should 
expect that our calculated superconducting critical tem- 
peratures may be somewhat overestimated compared to 
the experimental ones since we considered that all the 
neighboring tubes of a given one are metallic. For a more 
realistic description, we can consider a random distribu- 
tion of the tubes with different hclicities and diameters. 
Such complication is, actually, irrelevant for the physics 
of superconductivity in ropes of CNT since the nature of 
electronic transport is essentially sensitive to the trans- 
verse coupling between the tubes which depends on the 
intra-tube disorder [9]. 

The superconducting order is stabilized in the rope via 
Cooper pair tunneling between tubes and inside a single 
tube. We denote by Ji and J2 the Josephson coupling 
parameters across the rope, respectively, to the first and 
to the second neighboring tubes. We assume that the 
superconducting phase inside a tube is inhomogeneous 
with superconducting domains separated by metallic 
regions. This inhomogeneous structure, which may 
arise in the presence of impurities, is consistent with 
the absence of a bulk superconductivity in SWNT [21]. 
The superconducting domains along the tube (z axis) 
are coupled by Josephson tunneling parameterized by Jo- 

Regarding the strong superconducting fluctuations 
which extend on a large temperature range around the 
critical temperature T c , the mean field theory breaks 
down and one should expect clear deviation from the 
mean field critical temperature To. These fluctuations 
can be treated in the frame of the time dependent 
Ginzburg-Landau (TDGL) theory which has proven to 
be a reliable tool to study the critical transition region 
including superconducting fluctuations in different sys- 
tems such as high-T c [27] and low dimensional organic 
superconductors [28]. 

We start by writing the superconducting free energy F s 
of the rope compared to that of the normal state F norm : 



THE MODEL 

We consider a rope of identical SWNT arranged in a 
triangular lattice characterized by the basis (a, b) . For 
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where i and j denote the tube coordinates in the 
triangular basis (a, b) whereas n indicates the position 
of the superconducting domain along the tube direction 
z. ip n ij is the superconducting order parameter and () 
and (( )) correspond to the first and second neighboring 
tubes. The coefficients a and b are given by: a = a e 
and b = /j, n 2 elh 2 /2m 2 , where a = fr 2 /2m£ 2 , £ being 
the superconducting coherence length, and e = In (T/T ) 
while k = is the GL parameter. Here Ay and £y 
are respectively the London penetration depth and 
the coherence length in the (a, b) plane transverse to 
the tube direction. We take for simplicity £o = 
The Cooper pair is characterized by its electric charge 
eo = 2e and its effective mass m = 2m e where e is the 
unit charge and m e is the electron mass. 
The superconducting order is assumed to develop 
inside a tube over a thickness r 2 — v\ from the surface. 
The length of the superconducting domain is denoted Zqi • 

The Josephson parameters are written as: 



arp(— p), Ji = 



h 2 l e 
and Ji = - — —p2 ex P(~~k)> 
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where m* is the effective pair mass in the supercon- 
ducting domain whereas l e is the mean free path along 
the tube. L and D are, respectively the length and the 
diameter of the rope. We assume for simplicity that all 
the tubes have the same diameter. l\ (Z 2 ) denotes the 
intertube distance, from the tube surface, between first 
(second) neighboring tubes while l 02 is the distance be- 
tween superconducting domains inside a single tube. 

The natural question which arises concerns the origin 
of these Josephson coupling expressions. The major is- 
sue regards the exponential terms which lead to an en- 
hancement of the Josephson tunneling by increasing the 
amount of disorder, namely by decreasing the mean free 
path. 

This idea is based on previous studies dealing with 
Josephson coupled arrays of n-leg spin ladders [29] and 
correlated stripes in cuprate superconductors [30] which 
show clear evidence of the drastic effect of disorder on 
the superconducting state. Kivelson et cd. [30] have 
argued that the Josephson coupling between stripes is 
strongly enhanced by the transverse stripe fluctuations, 
which promotes the superconducting order. These fluc- 
tuations bring neighboring stripes close together leading 
to the enhancement of the mean value of the Josephson 
coupling. 



Orgad [29] has shown that such geometrical fluctua- 
tions in coupled ladder systems can reduce the suppres- 
sion of the superconducting correlations due to disorder, 
by increasing the Josephson tunneling between ladders. 
The dynamic of the ladders reduces the effective disor- 
der strength and make the superconducting pairing more 
robust against disorder. The interladder Josephson cou- 
pling is found to increase exponentially with the square of 
the fluctuation amplitude, which enhances the supercon- 
ducting transition temperature. Orgad [29] considered a 
Josephson tunneling amplitude depending on the inter- 
ladder distance as Jy ~ Joexp[— (s + ti, — Uj)/j], where 
Ui and Uj are the deviation of the i th and the ] th ladders 
from their static position, s is the mean distance of the 
ladder array and 7 is a characteristic constant [29]. 
The basic idea highlighted in Refs.[29, 30] is that the in- 
terplay between disorder and the dynamics of the stripes 
or the ladders is substantial for the stability of the su- 
perconducting order in cuprates and spin ladder super- 
conductors. 

Keeping this result in mind, let us now return to the 
rope of CNT. The latter can be described, as proposed 
by Ferrier et al. [9] by an array of ID atomic chains 
lying on a cylinder where each chain corresponds to a 
SWNT. The hopping processes along the chain are ran- 
domly distributed around a mean value t\\ with a square 
distribution 5tu . Such bond disorder along the chain may 
be induced by the dynamics of the tube as in the case of 
arrays of spin ladder or stripes. This leads to a competi- 
tion between the geometrical fluctuations of the SWNT 
and the local disorder inside the tubes. 

By analogy with Ref. [29] , the Josephson tunneling be- 
tween tubes can be written as J a exp[—d%j/'y], where 
dij is the separation distance between the \ th and the 
j*' 1 tubes. The exponential term expresses the Cooper 
pair tunneling probability which can be averaged over 
the tubes as (P±) = exp[—(d±)/ r y], where (d±) is an av- 
erage distance between the tubes. 

In diffusive superconductors, one should expect a de- 
pendence of the Josephson couplings on the mean free 
path since the superconducting coherence length is gov- 
erned by the disorder amount and reads as £ c = \J hv ^ 1 " , 
where A is the superconducting gap and vf is the Fermi 
velocity [31]. A key question raises at this point concern- 
ing the relationship between (dj_) and the intratubc mean 
free path l e , which we try to answer in the following. 

The plane transverse to the rope direction can be re- 
garded as a dirty two dimensional superconductor of a 
mesoscopic size where the disorder points, due to defects 
or impurities, are localized inside the tubes. In this plane, 
the tube sections form a sort of disordered clusters em- 
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bedded in a free disorder medium. The average distance 
between these clusters is controlled by the dynamic 
of the tube which is strongly dependent on the disorder 
amount inside the tubes. In the diffusive regime, the 
bond disorder due to the geometrical fluctuations of the 
tubes gives rise to an increasing intertube one particle 
hopping integral with increasing the site disorder ampli- 
tude originating from impurities and defects inside the 
tube [9]. This means that the intertube distance (d±) 
decreases with decreasing the intratube mean free path 
l e . (dj_) is then expected to have the same behavior as 
l e and may be expressed as a growing function of l e . We 
do not claim that the present model provides the exact 
form of this function. A more detailed analysis based on 
a microscopic study is needed. 

Since (d±), as l e , is a free parameter in our model, we 
set for simplicity (d±) = l e . This means that, in the 
diffusive regime, the mean free path inside the tube and 
across the rope are of the same order. This is justified as 
far as l e is smaller than the rope diameter D to keep the 
transverse one particle transport in the diffusive regime. 
Actually, this approximation does not affect the overall 
outcomes of our model but may yields to somewhat larger 
superconducting critical temperatures compared to the 
experimental ones. 

To characterize the electronic transport in disordered 
mesoscopic systems, one need to compare the size of the 
system, which is the rope diameter in this case, to a 
characteristic mean free path. Regarding its dependence 
on the intratube disorder amplitude, (d±) seems to be 
a good parameter to account for the transport regime 
across the rope. It comes out that (d±) and the rope 
diameter D, which depends on the tube number N, are 
the key parameters for the one particle transport and for 
the Cooper pair tunneling across the rope in the diffusive 
regime. The tunneling probability can then be written 
as (Pi) = exp[—(d±}/j] = exp[—l e /D], where the 7 con- 
stant, which accounts for the environment between the 
tubes, is replaced by rope diameter D. This is made 
possible since the tube environment is disorder free and 
depends only on the tube number included in the expres- 
sion of the rope diameter D. 

In the absence of site disorder and geometrical fluc- 
tuations, namely in a pure static rope, the Josephson 
couplings between respectively the first and the second 
neighboring tubes write as: 



Such couplings cannot describe the superconducting or- 
der in the rope since they are independent on the rope 
characteristics particularly the tube number. 
In the presence of disorder and geometrical fluctuations 
of the tubes, the Josephson parameters J\ and J2 given 
by Eq.3 should be changed to account for the average pair 
tunneling probability across the rope (Pi) = exp[— l e /D], 
which gives rise to the expressions introduced in Eq.2. 

Regarding the intratube Josephson tunneling Jo, one 
can define an average pair hopping probability along the 
tube (P||) = cxp[— l e /L] resulting from the geometrical 
fluctuations of the tube which yields to the expression 
given by Eq.2. 

It is worth to note that the Jo term is irrelevant for the 
stability of the superconducting phase as wc will show in 
the next. 

It comes out that the dynamics of the tubes in the 
rope mitigate the drastic effect of the local disorder on 
the superconducting order by enhancing the Josephson 
tunneling amplitudes between the tubes. The latter 
increase as a function of the effective disorder. This is 
reminiscent of the disorder-induced electronic transverse 
derealization in ropes of CNT proposed by Ferrier et al. 
[9]. We suggest that this derealization scenario holds 
for Cooper pair due to the tube dynamics as argued 
above. 

Let us now turn to the superconducting order param- 
eter whose critical dynamics satisfy the TDGL equation: 



r, 



dt 



OF 



+ (nij{r,t) 



(4) 



Here r o 1 = 7rft 3 /I6m£y k B T is the relaxation rate of 
the order parameter whereas £ n jj (r , t) are the Langevin 
forces describing the thcrmodynamical fluctuations and 
which obey the Gaussian white-noise law[27]: 

{Uj{r,t)Cn-i'A f ''*')) = 2T^k B TS(f- r')5(t - t') 

withr = (X + id,Y+jd,Z + nlo) and?"' = (X + i'd,Y + 
j'd, Z+n'lg), where d = li+do and Iq = ^01+^02, being 
the tube diameter. X, Y and Z are the coordinates of a 
point belonging to a superconducting domain of a SWNT 
of length Iqi, along the z direction, and of a thickness 

ri - n ■ 

By taking the derivative of the free energy (Eq.I) with 
respect to ^j*^, the TDGL equation becomes: 



Ji 



and Jo = 



2m*ll 
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(nij [r, t) = T — — h aip n!i j 
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I 
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where we adopted the Hartrcc approximation for the 
quartic term as in Rcf.27, which results in replacing the 
term b\ip ni j\ 2 ip ni3 by b(\ip n ij -\ 2 )tpnij ■ This approximation 
leads to a linear problem with a reduced temperature: 



a 



mj | 



which is determined sclf-consistently together with 
(IV'raijl 2 )- The superconducting critical temperature is 
defined as e(T = T c ) = 0[27]. 

To solve this equation, we introduce the Fourier trans- 
form of ipmj as 



where 



/•X2 pY~2 

ip(k,t) = V / dX dY 

nij Jx i Jy i 



dZ x 



( 6 ) ip nij {X + id,Y + jd, Z + nl , i)e( x+ "^ c (y+^)k y c (z+ni )k, ^ 



where X\ and X2 (Yi and I2) are the limiting values for 
X (Y) in a superconducting domain in the (a, b) plane. 
Taking the Fourier transform of Eq.5, we obtain 



"finij (r, t) 



d 3 k 



C(k, t) = I + 2m 7 + " + ^ _ cos(fc 2 Zo)) + 2Ji [3 - cos(dk x ) - cos(dk y ) - cos (d(k x - k y ))\ 

+ 2J 2 [3 — cos (d(2k x — k y )) — cos (d(k x + k y )) — cos (d(k x — 2k y )))} ijj(k,t), 

I 



(7) 



with a = a + b(\ip n ij\ } and the correlation relation sat- 
isfied by t) : 

(C(k,t)C(k',t')) = 2TQ 1 k B T{2Tif5{k - k')5{t - t 1 ) 

Equation 7 can be solved using the Green function 
method proposed by Puica and Lang [27] for lay- 
ered superconductors. We define the Green function 
R(k,t,k' z ,t') through the relation: 
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a\ = a + 



2m* 



d h 2 k 2 

r o 1 7T, + TT— T + 2J o (1 - cos(k z l )) + 01 
ot 2m* 



i?(fc, t, k' z ,t') = 5(k z - k' z )5{t - t'), 
where 



+ 2 J\ [3 — cos(dA; x ) — cos(dk y ) — cos (d(foc — k y ))] 



2J 2 [3 — cos (d(2k x — k y )) — cos (d(k x + k y )) — cos (d(k x — 2k y ))\ 

I 
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(9) 



We also introduce the Fourier transform of R(k,t,k' z ,t') 
with respect to time as: 

R{k,w,k' z ,t') = J dbRfat^J)^*-^ , (10) 

which can be deduced from Eq.8 as: 

R(k,u,k' z ,t') = S(k z - k' z ) x 



h 2 k 2 

-iwTq 1 + — -f + 2J (1 - cos(Mo)) + ai 
2m* 



,(H) 



tf>(k, t) solution of Eq.7 can be expressed in term of the 
Green function R(k,t, k' z ,t') as [27] 



ip(k,t) = I dt' I dk' g R(k,t,k' z ,l/)C(k x ,ky,k' g ,tf). 
Given Eq.ll, we obtain: 
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dr((k, t-r) <Ljc iut x 



h 2 kl 



-iujT 1 + — -4 + 2 J (1 - cos(Mo)) + ai 
2m* 



(12) 



with t — t — t' and the following correlation relation: 

(il>(k,t)i[>*(k',t)) a6(k-k') (13) 

To solve Eq.6, one need to derive (\ipnij\ 2 ) which, regard- 
ing Eq.13, can be simply written as: 



100 

N ( tube number) 



Wm 3 \ ) = 47rr - 1 fc B T / duj I d*k x 



2m 
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* + 2 J (1 - cos(k z l )) + oi 



(14) 



The critical temperature can now be deduced from 
Eq.6 by setting e(T = T c ) = 0, which yields to Eq.A.l 
given in the Appendix. In the next we discuss the nu- 
merical results. 



RESULTS AND DISCUSSION 

We have solved numerically Eq.A.l and the results arc 
depicted in Fig. 1 which shows the superconducting tran- 
sition temperature T c as a function of the number N of 
the tubes forming the rope. It is worth to note that iV 
is involved in the rope diameter D as D = s/N (da + e) 
where do and e are respectively the tube diameter and 
the intertube distance [13]. 

As shown in Fig.l, T c is strongly enhanced by increas- 
ing N but this enhancement is slowed down for N larger 
than 100 with a tendancy to saturation, which is remi- 
niscent of the experimental results [5, 26]. This behavior 
reflects the dimensionality of the superconducting phase 
appearing in the rope. By increasing N, the 3D char- 
acter of the superconducting state is enhanced and T c 
likewise. However, for a larger N (N ~ 200), the rope 
can be regarded as a 3D system and a further increase 
of N is irrelevant for the superconducting order, which 
explains the saturation behavior of T c at large N. 

A worth noting question concerns the interplay be- 
tween superconductivity and the ID character of a 
SWNT. Could superconductivity prevail over the low di- 
mensionality of such systems? This turns out to consider 
only the Jo term in our model. In such case, numerical 
calculations show that T c is at most of the order of 1 
mK, which explains the difficulty to observe an intrinsic 
superconductivity in SWNT as reported in Refs.9 and 
19. The superconducting phase can, actually, develop 
in ropes containing about one hundred metallic tubes as 



FIG. 1. Superconducting transition temperature as a func- 
tion of the number iV of tubes. The calculations are done 
in the one-particle delocalized regime and for A = 0.6/^m, 
£o = O.l^tm and L = 1.4/im [19]. A and £o are respectively 
the penetration depth and the coherence length in the super- 
conducting domain while L is the rope length. 



shown by earlier studies [16, 18]. The limiting tubes num- 
ber in our model is then N=13 if one include the first and 
second neighbors of a given tube. 

In Fig. 2, we give, for different tube numbers, the de- 
pendence of T c on the inverse of the mean free path which 
mimics the amount of local disorder inside the tube. Pe- 
culiarly, Fig. 2 shows that disorder promotes the super- 
conducting order as found experimentally [9, 19]. This 
behavior is due to the intertube disorder-induced delocal- 
ization of the Cooper pairs. Actually, the intertube pair 
delocalization is expected to develop in the electronic dif- 
fusive regime, where disorder can induce transverse hop- 
ping processes across the rope [9]. 

It is worth to note that the values of the critical tem- 
perature reported in Figs.l and 2 may be somewhat over- 
estimated since we have considered that all the tubes are 
metallic. In a more realistic model, one should take, on 
average, for each tube two neighboring metallic tubes 
since, in most cases, ^ of the tubes within a rope are 
metallic. 

By the way, one should emphasize the role of the 
Josephson tunneling J2 between second neighboring 
tubes on the superconducting order. Numerical results 
show that T c is reduced by 20% if J2 is neglected. Actu- 
ally, the second neighboring tubes should be involved in 
the tunneling processes since they are in the same range 
of reach as the first neighbors [20]. This is due to the 
geometry of the rope characterized by a tube diameter 
of 3 nm and an intertube distance of 0.35 nm. 

A worth stressing question regards the saturation be- 
havior of T c at large disorder amplitude in Fig. 2. This 
feature, which is due to the expression of the intertube 
couplings given by Eq.2, does not sound in agreement 
with the experimental data which rather show a collapse 
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FIG. 2. Superconducting transition temperature as a function 
of the inverse of the mean free path in a rope of SWNT for 
different tube numbers. The calculations are done in the one- 
particle delocalized regime and for the same data as in Fig.l 

of the superconducting phase at large enough amount 
of disorder[13, 26]. This discrepancy originates from 
the nature of the electronic transport regime. Our re- 
sults are derived within the delocalized diffusive regime 
characterized by a disorder induced transverse electronic 
hopping [9]. However, in the large disorder range, a lo- 
calized regime develops where the electrons are confined 
within individual tubes. The Josephson couplings given 
by Eq.2 are no more reliable since, in this case, the in- 
tratube disorder overcomes the geometrical fluctuations 
of the tubes, leading to the suppression of the intertube 
pair tunneling. The latter is expected to be strongly re- 
duced by the electron localization which can be roughly 
described by an expf— jf^ behavior for the intertube elec- 
tron hopping where L is the rope length and £ = 2Nl e 
is the localization length [9]. N and l e being the num- 
ber of metallic tubes and the mean free path inside the 
tube. As a consequence, one can assume the following 
Josephson couplings: 



which express the disorder induced Cooper pair localiza- 
tion as a result of the electronic localization. 

Fig. 3 shows the superconducting transition tempera- 
ture T c as a function of the inverse of the mean free path 
l e which is a measure of the disorder amplitude. The 
calculations are done using Eq.15. In this regime of lo- 
calization, T c is reduced by increasing disorder due to 
the suppression of the intertube tunneling. However, the 
tube number N acts, as in the delocalized regime, to 
the benefit of the superconducting phase. Increasing TV 
furthers the establishment of a 3D electronic transport 
regime by increasing the localization length £. The ef- 
fect of disorder is significantly important in ropes with a 




i , 1 , 1 , 1 , 1 , 1 

0,2 0,4 0,6 0,8 1 

Lfl (x 10" 3 ) 

FIG. 3. Superconducting transition temperature as a function 
of the inverse of the mean free path in ropes of N=30 and 
N=70 tubes. T c is calculated in the localization regime where 
the Josephson couplings are given by Eq.15. The used data 
are the same as in Figs.l and 2. 

small tubes number where the ID character prevails over 
the formation of a 3D superconducting order. 

The superconducting behaviors in the delocalized and 
localized regimes (Figs. 2 and 3), are reminiscent of those 
obtained in a 2D array of stripes [30]. In such systems, 
the superconducting transition temperature is found to 
increase with the transverse stripe fluctuations up to a 
critical value above which it drops. This happens when 
the system undergoes a phase transition to an isotropic 
state where the stripe structure is lost. 

A tough question raised from Figs. 2 and 3 concerns 
the extension of the delocalized regime. At which dis- 
order amplitude the dynamic of the tubes is frozen and 
the intertube Josephson tunnclings start to collapse? A 
rough estimation may be deduced from the experimen- 
tal results of Kasumov et al. [13] showing that the key 
parameter governing the disorder in a suspended rope is 
the ratio *f- where L and £ c are respectively the rope and 
the coherence lengths. The latter depends on the mean 

free path l e as discussed above £ c = J %v ^ e [31]. 

In Fig. 4, we have depicted the behavior of supercon- 
ducting transition temperature in the localized and de- 
localized regimes for a rope of N=70 tubes based on the 
results shown in Figs. 2 and 3. According to Fig. 4, the 
suppression of the superconducting order starts at a crit- 
ical value y~ =0.2. The smaller the tube number, the 
greater l ec , the frailer the superconducting order. 

The comparison of the numerical values of l ec with 
the experimental results of Ref.9 is not obvious. More 
data are needed to accurately determine the critical dis- 
order amplitude at which the superconducting transi- 
tion temperature reaches its maximum before decreas- 
ing. Nevertheless, one can compare the extent of the dis- 
order regime over which the superconducting order de- 
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0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 
L/f (xlO" 3 ) 



FIG. 4. Superconducting transition temperature as a function 
of the inverse of the mean free path in ropes of N=70 tubes. 
Regions (I) and (II) denote, respectively, the delocalized and 
localized regimes. The calculations are done with the same 
data as in Fig.l 



on the stability of the superconducting order in ropes of 
CNT. We found that an increase of the dimensionality of 
the rope, which is achieved by increasing the tube number 
N, promotes the establishment of a 3D superconducting 
phase with an increasing superconducting critical tem- 
perature T c . However, for large N values, T c tends to sat- 
uration indicating the formation of a well defined 3D su- 
perconducting order. The main result of our work regards 
the disorder induced superconductivity in the rope which 
originates from the dynamics of the tubes. The latter en- 
hance the intertube Josephson tunnelings which mitigate 
the suppression of the superconducting phase by disor- 
der. However, for larger disorder amplitude, electronic 
localization prevails against intertube hopping leading to 
the suppression of superconductivity as found in other 
superconducting materials. 
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APPENDIX: SUPERCONDUCTING CRITICAL 
TEMPERATURE 

By setting e(T = T c ) = in Eq.6 we obtain the follow- 
ing equation giving the superconducting critical temper- 
ature: 



CONCLUDING REMARKS 

In summary using TDGL theory, we probed the role of 
the effective dimensionality and the amount of disorder 



J 



, Tc r 

ln 77T + 9 T / sm 
Jo Jo 



2tt r c r W c 

2, 



dip k z dk \ dWf(k,e,(p,W) = 0. 
Jo Jo 



(A.l) 



The /(fc, 9, ip, W) is given by: 



f(k, 9,tp, W) = { W 2 + ^ + 2^ [1 - cos(£:Z cos 9)] 
2m* ciq 
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2— [3 — cos(kd sin 9 cos ip) — cos(kdsin6sm<p) — cos (kdamd(cos(p — sin </?))] 

2— [3 — cos (fco?sin6'(2 cos ip — sin <p)) — cos (/cdsin^fcos y> + sin ip)) — cos (fee? sin 6*(cos (p — 2 sin p))](A.2) 
a 



r 



We introduced the dimensionless variable W 

,2„2 c4 



tHlj 



8k B T 

and we set g = 87r/io^B K2e o£|[ > with A = 0.6/im and 

=e = 0.1^ [i9]. 

We have adopted the no-cutoff limit for c and W c (c — > 
oo, W c — ► oo), which means that all types of supercon- 
ducting fluctuations, even with short wave lengths, arc 
considered [27]. 
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